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CROSS-REFERENCE TO RELATED APPLICATIONS 
This application Is based upon and claims the 
benefit of priority from the prior Japanese Patent 
Application No. 2003-308919, filed on September 1, 
2003, the entire contents of which are Incorporated 
herein by reference. 

BACKGROUND OF THE INVENTION 
[Field of the Invention] 

The present Invention relates to a method and 
equipment for simulation which correct flare (local 
flare) locally occurring In the photo process of 
semiconductor manufacturing. 

[Description of the Related Art] 

Pattern formation technology In manufacturing a 

semiconductor device adopts optical lithography. The 
optical lithography refers to a method In which a 
desired pattern formed In a reticle as a photo mask Is 
transferred to a sample substrate through a reduction 
optical system. Photosensitive resin called a resist 
Is formed on the sample substrate. With the use of 
difference In the melting speed between exposed areas 
and unexposed areas by developing solution, a latent 
Image of the transferred pattern Is etched to desirably 
process material. 

Exposure technology adopts optical proximity 
correction (OPC), In which the effect of traces In 
proximity to each other Is calculated and corrected to 
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transfer the fine pattern with high precision. In 
addition to optimizing the wavelength of exposure light 
and the structure of the reticle. A transferred Image 
of the pattern Is simulated in consideration of the 
Illumination condition (NA and Sigma) and the exposure 
condition (the material of the resist and exposure 
wavelength) and the like of an exposure device, and the 
effect of optical proximity and correction values are 
calculated or determined by experiment, to correct the 
dimensions of the reticle. 

Generally, when the optical proximity correction 
is carried out, an optical image of the layout of the 
photo mask (mask layout) is calculated by a simulator, 
and then the mask layout is so modified that the 
calculated optical image takes a predetermined shape. 
It is necessary as a precondition that the optical 
image outputted from the simulator precisely 
corresponds with actual one. The unevenness in the 
material of a lens of the exposure device and the 
specific characteristics thereof, however, cause 
scattering in the lens, so that the effect of 
scattering also has to be considered in the simulation. 
The effect of scattering, called local flare, is caused 
by the specialty of lens material which depends on the 
wavelength of exposure light (short wavelength 
represented by 193nm) and the like. The effect of 
scattering becomes a main cause of unexpected 
deformation in the shape of a transferred pattern and 
in the width of traces. In an image surface, it is 
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possible to approximate the local flare at even light 
Intensity In a background which depends on the 
Intensity transmlttance of the photo mask. 

Patent document 1 : Japanese Patent Application 
Laid-open No. 2001-272766. 

In a conventional simulation method, when the 
dimensions of traces are enough larger than an 
exposure wavelength. It Is possible to estimate the 
amount of local flare with sufficient accuracy. When 
there are traces the dimensions of which are equal to 
or smaller than the exposure wavelength, however, 
there Is a problem that the amount of local flare 
cannot be obtained with high accuracy. As the 
pattern of a semiconductor device becomes finer and 
finer in recent years. It Is desirable to develop a 
simulation method which estimates the amount of local 
flare with high accuracy even if a trace is fine. 

SUMMARY OF THE INVENTION 
Considering the foregoing problem, an object of 

the present invention is to provide a method and 

equipment for simulation which estimate the amount of 

local flare with extremely high accuracy and precisely 

predict the occurrence of the local flare, not only 

when the dimensions of traces are larger than an 

exposure wavelength, but when the traces equal to or. 

smaller than the exposure wavelength are mixed in a 

pattern . 
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A method for simulation according to the present 
Invention, which simulates an amount of occurrence of 
local flare occurring In an exposure process In 
manufacturing a semiconductor device, comprises the 
steps of dividing a layout of a photo mask Into a 
*v plurality of areas, calculating an average value of 

light Intensity In each of the areas, and estimating 
the amount of occurrence of local flare In each of 
the areas on the basis of each of the average values. 

Equipment for simulation according to the present 
Invention, which simulates an amount of occurrence of 
local flare occurring In an exposure process In 
manufacturing a semiconductor device, comprises 
division means for dividing a layout of a photo mask 
Into a plurality of areas, and average light 
Intensity value calculation means for calculating an 
average value of light Intensity In each of the areas. 
In order to estimate the amount of occurrence of 
local flare In each of the areas on the basis of each 
of the average values . 

According to the method and equipment for 
simulation of the present Invention, even If traces 
equal to or smaller than an exposure wavelength are 
mixed In a pattern. It Is possible to estimate the 
amount of local flare with extremely high accuracy, 
and precisely predict the occurrence of the local 
flare. The dimensions of traces are easily and 
precisely corrected based on simulation results, so 
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that It Is possible to actualize a highly reliable 
semiconductor device . 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. lA Is a schematic plan view showing the mask 
layout of a photo mask In which the amount of local 
flare Is evaluated In a conventional method, and Figs. 
IB and IC are explanatory views showing an example of 
Intensity transmit tance and an example of evaluation 
values ; 

Fig. 2 Is a block diagram showing the schematic 
structure of simulation equipment according to the 
present Invention ; 

Fig. 3 Is a flow chart showing a simulation 
method according to the present Invention; 

Fig. 4A Is a schematic plan view showing the mask 
layout of a photo mask In which the amount of local 
flare Is evaluated, and Figs. 4B and 4C are 
explanatory views showing an example of light 
Intensity distribution and an example of evaluation 
values ; 

Fig. 5 Is a schematic diagram which explains a 
method of determining the weighting factor of 
diffracted light; 

Fig. 6 Is a schematic diagram which explains a 
method of determining the weighting factor of 
diffracted light, when a ring-shaped lamp Is used; 

Fig. 7 Is a block diagram showing the Internal 
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structure of a personal user terminal; 

Fig. 8 Is a schematic plan view showing the mask 
layout of a photo mask according to a first 
embodiment ; 

Fig. 9 Is a graph showing the relation between a 
line pitch and a light Intensity value In the photo 
mask ; 

Fig. 10 Is a schematic plan view showing the mask 
layout of a photo mask according to a second 
embodiment ; and 

Fig. 11 Is a graph showing the relation between a 
position In the photo mask and a light Intensity 
value . 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
-ESSENTIAL FEATURES OF THE PRESENT INVENTION- 
Flrst, the essential features of the present 
Invention will be described. 

The Inventor was dedicated to studying a method 
for estimating the amount of local flare with high 
accuracy, even If traces equal to or smaller than an 
exposure wavelength are mixed in a pattern, and has 
found that the method is built on the following 
theory . 

Fig. lA Is a schematic plan view showing the mask 
layout of a photo mask in which the amount of local 
flare is evaluated in a conventional method, and Figs. 
IB and IC are explanatory views showing an example of 
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Intensity transmit tance and an example of evaluation 
values . 

In the conventional method, as shown In Fig. lA, 
the mask layout Is divided Into a plurality of 
rectangular areas. The transmlttance (open ratio) of 
each area calculated and smoothed Is estimated as the 
amount of local flare. In other words, when an area 
(l,j) Is decided as shown In Fig. IB, the predicted 
amount of local flare F(l,j) In the area (l,j) Is 
obtained by the following equation (1). Fig. IC 
shows the calculation results of each area. 

,'=-2;'=-2 

Wherein, T{1,3) Is the average value of 
transmlttance In the area (1,3), and E(l',j') Is a 
weighting factor which represents the degree of 
effect of the area. For convenience of explanation, 
only the equation for the area (1,3) Is described 
above, and smoothing processing Is carried out with 
respect to two areas adjoining forward and backward. 

In the foregoing conventional method, there Is a 
case where simulation result does not match with the 
actual amount of local flare, because the amount of 
local flare Is estimated on the basis of the 
Intensity transmlttance of the mask layout before 
passing through a lens of an exposure device. In 
other words, light passing through a fine mask 



pattern in the photo mask Is diffracted, and light 
diffracted with a large angle does not pass through 
the lens. 

Considering that only light passing through the 
lens of the exposure device, out of light passing 
through the mask pattern. Is provided for exposure, 
the Inventor has come up with an Idea of using the 
average value of light Intensity distribution In an 
Image surface on which the light forms an Image of 
the mask pattern. Instead of the average value of the 
Intensity transmit tance of the mask pattern used In 
the conventional method In order to grasp only the 
effect of the above-mentioned light accurately. 

-CONCRETE EMBODIMENTS OF THE PRESENT INVENTION- 

Flg. 2 Is a block diagram showing the schematic 
structure of simulation equipment according to the 
present Invention, and Fig. 3 Is a flow chart of a 
simulation method according to the present Invention. 

The simulation equipment comprises division means 
1 for dividing the layout of a photo mask (mask 
layout) Into a plurality of areas, average light 
Intensity value calculation means 2 for calculating 
the average value of light Intensity In each area, 
smoothing means 3 for subjecting the calculated 
average value to smoothing processing, and 
multiplication means 4 for multiplying the smoothed 
average value by a predetermined multiplier. To 
evaluate the amount of local flare with the use of 
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the simulation equipment, first, as shown In Fig. 4A, 
the division means 1 divides the mask layout Into a 
plurality of areas (step 1). Then, the average light 
Intensity value calculation means 2 calculates the 
average value of light Intensity In each area (step 
2), and the smoothing means 3 smoothes the calculated 
average value of light Intensity (step 3). Then, the 
multiplication means 4 multiplies the smoothed 
average value of light Intensity by a predetermined 
multiplier (step 4), and each obtained value Is 
evaluated as the amount of occurrence of local flare 
In each area (step 5). 

To be more specific, as shown In Fig. 4B, the 
average value of light intensity 1(1, j) is used 
Instead of T(l,3) in Fig. 1. The amount of local 
flare F(l,j) in an area (1,3) to be predicted is 
calculated by the following equation (2). Fig. 4C 
shows the calculation result of each area. 

^0»;)=ZS^0'.7')x/(/+/',y+y') ~. (2) 

i'=-2;'=-2 

wherein, 1(1,3) is the average value of light 
intensity in the area (1,3), and E(i',3') is a 
weighting factor which represents the degree of 
effect of the area. For convenience of explanation, 
only the equation for the area (1,3) is described 
above, and the smoothing processing is carried out 
with respect to two areas ad3olnlng forward and 
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backward . 

According to this method, the amount of local 
flare Is calculated on the basis of the average value 
of light Intensity of an optical Image on a wafer 
after passing through a lens of an exposure device. 
By doing so, since light not passing through the lens 
Is not Included In the calculation of local flare, 
the amount of local flare Is calculated In a 
situation closer to actual one, so that predictive 
accuracy Is Improved. 

To carry out the foregoing calculation, it Is 
necessary to obtain the light Intensity of the mask 
layout. It is unrealistic, however, to calculate the 
light intensity of the whole mask layout by use of a 
conventional light intensity simulation method, for 
the reason of processing time. Thus, noting that 
what is necessary is the average value of light 
intensity, diffracted light is calculated by a 
Fourier transformed image in each area of the mask 
layout, and the average value of light intensity is 
obtained by multiplying the light intensity of the 
diffracted light passing through a projection lens by 
a predetermined multiplier. In the concrete, the 
average value of light intensity (I) is calculated 
with the use of the following equation (3). 

I=f,F,S,S: ... (3) 

k-l 
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wherein, "Fjc" Is a weighting factor of diffracted 
light (a real number), and "S^" Is the amplitude of 
the diffracted light (a complex number: an asterisk 
(*) represents a complex conjugate number. 
representing the number of the diffracted light 
corresponds to the combination of (n,m) one to one. 
("n" and "m" are the orders of diffracted light In X 
and Y directions, respectively.) 

The average value of light Intensity Is obtained 
area by area (1,J). On the assumption that the area 
(l,j) Is periodic, Fourier transforming the area 
(l,j) makes It possible to obtain the diffracted 
light . When the area Is enough larger than an 
exposure wavelength (approximately five times), there 
Is no problem with the assumption that the area Is 
periodic. The weighting factor "Fjc" of the 
dlffractlve light Includes Information whether or not 
the dlffractlve light passes through the lens. The 
weighting factor "F^" of the dlffractlve light Is 
obtained as shown In Fig. 5. 

Referring to Fig. 5, "C" Is the center of the 
lens In the exposure device, "Sk" Is the position of 
diffracted light, "NA" Is the numerical aperture of 
the lens, and "a" Is the radius of light source with 
respect to "NA" (an Interference factor). The area 
"Ale" shared between a circle "C and a circle "Sk" Is 
In proportion to the weighting factor of the 
diffracted light "Sk". The area "Ak" Is also In 



proportion to a probability that the diffracted light 
passes through the lens. "Ak" is divided by "a^Tc" and 
normalized to obtain the weighting factor "Fk" . The 
definitive weighting factor is represented by the 
following equation (4). 

F,=Aj{cT'7t) ... (4) 

When a ring-shaped lamp is used, the weighting 
factor is obtained as shown in Fig. 6, and is 
represented by the following equation (5). 

Pk = A /(o-2^^ - aiTt) ... ( 5 ) 

The circular or ring-shaped light source is used 

in the foregoing description, but the essence is to 

calculate the area shared between the shape of a 

light source and a lens circle. It is apparent that 

this method is applicable in a lighting situation 

except for above. 

The simulation value of the amount of local flare, 

obtained as described above, is used for optical 

proximity correction (OPC), or, as described later, 

is added to a light intensity value by optical 

simulation in order to provide it for the simulation 

of an optical image. 

According to this embodiment, as described above, 

it is possible to actualize simulation by which the 
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amount of local flare Is estimated with extremely high 
accuracy and the occurrence of the local flare Is 
precisely predicted, not only when the dimensions of 
traces are larger than the exposure wavelength, but 
when the traces equal to or smaller than the exposure 
wavelength are mixed In the mask pattern* 

By operating a program product stored in a RAM, 
ROM or the like of a computer, it is possible to 
actualize each means (means 1 to 4 of Fig. 2 and the 
like) for constituting the foregoing simulation 
equipment, and each step (the steps 1 to 5 of Fig. 3 
and the like) in the simulation method. The present 
invention includes that program product, and a 
computer-readable recording medium in which the 
program product is stored. 

To be more specific, the program product is 
provided for the computer via the recording medium 
such as, for example, a CD-ROM or various 
transmission media. As the recording medium for 
recording the program product, a flexible disk, a 
hard disk, a magnetic tape, a magneto-optical disk, a 
nonvolatile memory card or the like is available in 
addition to the CD-ROM. As the transmission medium 
of the program product, on the other hand, a computer 
network (WAN such as LAN, the Internet and the like, 
a wireless communication network and the like) or a 
communication medium in a system (a wire circuit such 
as an optical fiber and the like, a wireless circuit 
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and the like) Is available. 

Not only when the function of the foregoing 
embodiment Is actualized by running the program 
product provided for the computer, but when the 
program product actualizes the function of the 
foregoing embodiment In cooperation with an OS 
(operating system), other application software or the 
like which the program product runs on the computer, 
or when the whole or a part of processing of the 
provided program product Is carried out by a function 
enhancement board or unit of the computer to 
actualize the function of the foregoing embodiment, 
the present Invention Includes that program product. 

Fig. 7, for Instance, Is a block diagram showing 
the Internal structure of a general personal user 
terminal. Referring to Fig. 7, the reference numeral 
1200 denotes a personal computer (PC). The PC 1200 
having a CPU 1201 runs device control software stored 
In a ROM 1202 or a hard disk (HD) 1211 or provided by 
a flexible disk drive (FD) 1212, In order to control 
each device connected to a system bus 1204 as a whole. 

-EMBODIMENTS- 

Hereafter, some embodiments In which the amount 
of local flare Is actually estimated with the use of 
the method according to the present Invention will be 
described . 

( First Embodiment ) 

In the first embodiment, the average value 



% 

(relative value) of light intensity of a photo mask 
having a layout shown in Fig. 8 is obtained with the 
use of the method according to the present invention. 
The photo mask has a line-and- space pattern in which 
the width of an opening 11 is equal to that of a 
light shielding portion 12. It is assumed that the 
length "L" of the opening 11 is much longer than the 
width "W" thereof. 

Fig. 9 is a graph of characteristics showing the 
relation between a line pitch in the photo mask and 
the value of light intensity. 

Fig. 9 shows the average values of light 
intensity calculated by a method of the present 
invention, in relation to line pitches "P". The 
average values of intensity transmit tance according 
to a conventional method are also put therein for 
comparison purpose. In this embodiment, the average 
value is calculated on the condition that "NA" is 0.7, 
an exposure wavelength is 193nm, and an interference 
factor is 0.8. 

According to the conventional method, since the 
width of the opening 11 of the photo mask is equal to 
that of the light shielding portion 12, the average 
value of intensity transmit tance becomes constant 
(0.5) irrespective of the line pitch. According to 
this method, on the other hand, diffracted light 
passing through a lens decreases as the line pitch 
becomes smaller, so that the average value of light 
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Intensity decreases as a result. This means that 
this method reflects the variation In the amount of 
local flare depending on the line pitch, and produces 
realistic values. 

( Second Embodiment ) 

In the second embodiment, the amount of local 
flare calculated by the method of the present 
Invention Is added to an optical Image by optical 
simulation, to obtain an optical Image In 
consideration of local flare. In this embodiment. It 
Is assumed that the ratio of the amount of local 
flare to the Intensity of the optical Image Is 5%. 

Fig. 10 shows the shape of a photo mask used In 
the calculation of this method. The photo mask 
Identical to that of Fig. 7 has a line - and- space 
pattern In which the width of an opening 11 Is equal 
to that of a light shielding portion 12, and has a 
line pitch of 0.5|j.m. It Is assumed that the length 
"L" of the opening 11 Is much longer than the width 
"W" thereof. 

Fig. 11 Is a graph of characteristics showing the 
relation between the position In the photo mask and 
the value of light Intensity. The position In Fig. 
11 corresponds to the position In the photo mask 
along the broken line I-I of Fig. 10. 

Fig. 11 shows the average value of light 
Intensity calculated by the method of the present 
Invention and subjected to smoothing processing, the 
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optical Image by optical simulation, and the result 
of adding the average value of light Intensity to the 
optical Image, namely adding the amount of local 
flare of 5% to the optical Image. As described above, 
the average value of light Intensity calculated by 
this method contributes to the accurate simulation of 
an optical Image In consideration of local flare. 



